Conformal-thin-sandwich initial data for a 
single boosted or spinning black hole puncture 



Pablo Laguna 
Center for Gravitational Physics and Geometry, 
Center for Gravitational Wave Physics, 
Department of Astronomy and Astrophysics, 
Department of Physic, 
Penn State University, University Park, PA 16802 



cn 
o 
o 

(N 
-t— > 
O 

O 



Sequences of initial-data sets representing binary black holes in quasi-circular orbits have been 
used to calculate what may be interpreted as the innermost stable circular orbit. These sequences 
have been computed with two approaches. One method is based on the traditional conformal- 
transverse-traceless decomposition and locates quasi-circular orbits from the turning points in an 
effective potential. The second method uses a conformal-thin-sandwich decomposition and deter- 
mines quasi-circular orbits by requiring the existence of an approximate helical Killing vector. Al- 
though the parameters defining the innermost stable circular orbit obtained from these two methods 
differ significantly, both approaches yield approximately the same initial data, as the separation of 
the binary system increases. To help understanding this agreement between data sets, we consider 
the case of initial data representing a single boosted or spinning black hole puncture of the Bowen- 
York type and show that the conformal-transverse-traceless and conformal-thin-sandwich methods 
yield identical data, both satisfying the conditions for the existence of an approximate Killing vector. 

PACS numbers: 04.30+x 
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I. INTRODUCTION 

Binary systems involving compact objects such as 
black holes and neutron stars have been the main source 
of attention of numerical relativists in the past two 
decades. The motivation for studying these systems 
has been enhanced by the urgency of producing results 
that could help observational efforts like LIGO, GEO600, 
TAMA300 and VIRGO Q. Unfortunately the problem 
of producing complete and sufficiently accurate numeri- 
cal simulations of orbits and mergers of compact binaries 
is still not feasible. This impediment does not imply that 
the results produced so far are devoid of any astrophysi- 
cal relevance. Simulations of a few orbits of neutron star 
binaries have been obtained 0, Q • Grazing collisions and 
wave forms produced from the final plunge of a binary 
black hole system have been also computed Q. And of 
direct relevance to the work presented in this paper is the 
study of the behavior of binary orbits using sequences of 
initial data configurations. 

Astrophysically realistic initial data sets representing 
black hole binaries are crucial. Without them, the pre- 
dictive power of gravitational wave source simulations, 
with relevance to data analysis efforts, is seriously com- 
promised. Constructing initial data sets involves a se- 
ries of choices such as conformal transformations, tensor 
decompositions, topologies, boundary conditions, freely 
specifiable data to name a few. These choices have a 
profound effect on the physical properties of the data. 

Two general approaches for constructing initial data 
are currently under consideration. The oldest approach 
has its roots in the work by Lichnerowicz and York 
This approach is based on conformal transformations 
and transverse-traceless (CTT) decompositions. The sec- 



ond and most recent approach, currently receiving con- 
siderable attention, attempts to establish a more di- 
rect "space-time" connection. It borrows some of the 
ideas from the CTT method, namely conformal trans- 
formations and traceless decompositions, but in addi- 
tion provides a natural framework for setting up quasi- 
cquilibrium configurations of binary systems. Under this 
approach, initial data can be viewed as derived from a 
"thick slice or thin sandwich" section of the space-time. 
The method was pioneered by Mathews and Wilson [||; 
recently used by Grandclement, Gourgoulhon and Bonaz- 
zola ; and formally characterized by York || . Because 
the approach also involves conformal transformations, it 
is known as the conformal thin-sandwich (CTS) method. 

Initial data sets representing binaries in quasi-circular 
orbits will be needed as a starting point of any evolution 
simulation aimed at producing gravitational waveforms 
of astrophysical black hole binaries. Cook introduced 
the first method to locate what could be identified as 
black hole binaries in quasi-circular orbits. The method, 
called the effective potential method, is based on the po- 
tential 



E h = E — 2M . 



(1.1) 



where Eh represents the binding energy of the system, 
E the total energy and M the irreducible mass of each 
individual black hole. The total energy E is a well de- 
fined quantity, namely the total ADM mass ^(j- How- 
ever, defining individual black hole masses for interact- 
ing black holes is not a completely well defined concept 
[llj . Nonetheless, a reasonable approximation is to as- 
sume that the irreducible mass is given by the proper 
area of the apparent horizon. Quasi-circular orbits are 
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then found from the condition 
dE b 



dl 



M,J 



= o, 



(1.2) 



where I is the proper separation of the horizons and J 
the total angular momentum of the binary. Similarly, 
the angular velocity £1 is found from 



dE b 
dJ 



MA 



(1.3) 



Given this approach to identifying quasi-circular orbits, 
the innermost stable circular orbit (ISCO) is found by 
looking at the end point of a sequence of quasi-circular 
orbits. 

Baumgarte 0] and later Baker recomputed, us- 
ing Cook's effective potential method, quasi-circular or- 
bit sequences. One of the motivations for this work was 
to investigate whether Cook's results were sensitive to 
how the black hole singularities are handled. Instead of 
using a topology requiring conformal-imaging as used by 
Cook |[g, Baumgarte and Baker handled the black hole 
singularities via the puncture method [14( . Both studies 
produced results in close agreement with those by Cook. 
In the remainder of this paper, quasi-circular orbit data 
sets computed with the CTT method will be referred as 
Cook-Baumgarte-Baker (CBB) data sets. 

As mentioned before, the CTS approach has also 
been used by Gourgoulhon, Grandclement and Bonaz- 
zola (GGB) 0, llj| to construct quasi-circular orbits of 
black hole binary systems. The essence of this work is 
to identify quasi-circular orbits based on the assump- 
tion that the space-time possesses an approximate heli- 
cal Killing vector. One important aspect of the GGB se- 
quences is that they yield an ISCO identification in better 
agreement with post-Newtonian results jl6j | when com- 
pared with CBB sequences. Specifically, the normalized 
binding energy e = E/2 M — 1 at the ISCO from post- 
Newtonian estimates is e w —1.67%. The GGB study 
yields e i=s —1.73%. On the other hand, the CBB estimate 
is significantly more strongly bound, e ~ —2.3%. These 
differences show more dramaticaly when comparing the 
ISCO orbital period. Post-Newtonian calculations yield 
T/2M w 71.2 and in the GGB work T/2M « 62.8. In 
contrast, the CBB data estimate T/2M ps 37. 

Reasons for these discrepancies are not completely un- 
derstood. One suggestion was that these differences were 
due to the criteria for identifying quasi-circular orbits. 
Instead of the method used by CBB described above, 
GGB determined the co-rotating frame defining quasi- 
circular orbits from the requirement that the ADM and 
Komar 01 masses are equal. The motivation for this 
choice was that the Komar and ADM masses only agree 
for stationary space-times. Therefore, the data com- 
puted this way have a better chance of representing a 
quasi-stationary situation. Skoge and Baumgarte [l|| 
designed a simple example involving thin-spherical shells 
of collision-less particles in circular orbits which shows 
that the effective potential method and the equal Komar 



and ADM masses method for identifying circular orbits 
are equivalent. More recently, Tichy and Briigmann [l9l | 
computed sequences of CTS binary black hole punctures 
that yield a similar conclusion. These studies provide 
strong support to the view that the differences between 
CBB and GGB data sets are not due to the method for 
identifying quasi-circular orbits but somewhere else. 

Another possible explanation for the ISCO discrepancy 
that has been mentioned is once again the approach to 
representing the black hole singularities. The GGB ini- 
tial data sequences are based on a topolo gy requiring 
conformal-imaging 0- The work in Refs. |l9l |20| was 
mostly aimed at testing this by extending the GGB re- 
sults to the case of black holes represented by punctures. 
Unfortunately, the data obtained in this work only ap- 
proximately satisfied one of the conditions used by GGB 
to define quasi-stationarity. Recently, Hannam and col- 
laborators j2]| pointed out the serious difficulties in con- 
structing CTS binary data in quasi-equilibrium based on 
punctures instead of conformal-imaging. 

There is however on important aspect of the CBB and 
GGB sequences that cannot be ignored. In spite of their 
differences for determining the ISCO, the CBB and GGB 
sequences show a close agreement when the separation of 
the binary systems increases. Therefore, at large separa- 
tions, it should be possible to translate CBB into GGB 
data and vice- versa, independently of the procedure used 
to represent the black holes. The work in this paper is 
motivated by this observation. Since for large binary 
separations the black holes are approximately isolated, 
we focus attention on single boosted or spinning black 
hole punctures. We demonstrate that it is possible to 
translate the CTT data into a CTS form while preserv- 
ing the conditions used by GGB, namely the existence 
of an approximate Killing vector, thus circumventing in 
these cases the problems pointed out in Ref. [2l| . 

The paper is organized as follows. In Sec. [H] the con- 
formal decomposition of the Einstein constraints is pre- 
sented up to the point where the CTT and CTS methods 
agree. Brief reviews of the CTT and CTS methods are 
presented in Sees III Al and III Bl respectively. The con- 
ditions necessary for translating between CTT and CTS 
data are derived in Sec. IIIII The main results regarding 
boosted or spinning punctures are presented in Sec. IIVI 
Conclusions are given in Sec. 



II. DECOMPOSITION OF THE CONSTRAINTS 



Under the standard 3+1 decomposition of the Einstein 
equations, the 4-dimensional line element is written as 

ds 2 = -a 2 dt 2 + g tj (dx l + ft dt) {dx 3 + ft dt) , (2.1) 
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with cjij the intrinsic 3-metric to S t (t = constant hyper- 
surfaces), a the lapse function and /3 1 the shift vector. 1 
The extrinsic curvature of S t is defined by 

k lj = --C N g lj , (2.2) 

with Cm the Lie derivative along the unit time- like nor- 
mal to St- In terms of these variables, the vacuum 
Einstein equations consist 01 °f a set °f evolution equa- 
tions: 

c?« ffij = -2aK i:j (2.3) 
d Kij = —ViVja + aRij 

+ a{KK lJ -2K lk K k 3 ), (2.4) 

and a set of constraint equations: 

R + K 2 - k i:i K l: > = (2.5) 
VjiT^' - V l J^ = . (2.6) 

Above, Vi and Rij denote respectively covariant dif- 
ferentiation and Ricci tensor associated with g^. Also 
K = g^Kij, R = g^Rij and o = d t — Cp. 

Constructing initial data involves specifying (gij, K^) 
on a given S t subject to the four constraint equations 
(|2.5|) and 1)2. 6fl . In other words, out of the twelve pieces 
in (jjij, Kij), only eight are freely specifiable. The re- 
maining four are fixed by these constraints. Thus, setting 
up the initial data problem in general relativity becomes 
the task of singling out in (g^ , Kij ) those quantities that 
will be fixed by the constraints. This is accomplished 
via conformal transformations and tensorial decomposi- 
tions. The CTT and CTS approaches are examples of 
this methodology. 

Both the CTT and CTS methods involve applying a 
conformal transformation of the 3-metric and a decom- 
position of the extrinsic curvature of the form: 

9ij = <t> 4 9ij (2-7) 
K^ = A^ + ^g^K, (2.8) 

such that A % i = 0. Notice that the choice K — K has 
been made. In terms of these new variables, the con- 
straint equations ((2.5(1 and 1)2.6(1 become: 

cj)- A R - 8ct>- 5 V 2 cl>+~K 2 - AijA ij =0 (2.9) 

O 

VjA ij + 10A lj Vj\ncj)~?-<j)- 4 V t K = 0,(2.10) 

o 

where Vj and R denote respectively covariant differen- 
tiation and Ricci scalar associated with the conformal 



1 Latin indices denote 3-dimensional spatial indices and Greek in- 
dices 4-dimensional space-time indices. 



3-metric g^. In deriving Eq. 1(2.10(1 . we have used the 
following identity: 

Vj f ij = <t> n [VjT 13 + (10 + n) T ij Vj In <j>\ , (2.11) 

where T 4 -? = (f) n , with n an exponent. Similarly, the 
evolution equations 12.3(1 and 1(2.4(1 can be rewritten as: 

d t <t> = ft Vi + i (VijS* - a K) (2.12) 
dtgij = (L/3)y -2a0- 4 iy (2.13) 
d K = -V 2 a + a [A^ A l} + ^K 2 ^j (2.14) 

d Ay = [~V 1 V J a + aR, j ] TF 

- a(2A\ k A k j + ^KA^ (2.15) 

where the superscript TF denotes the trace-free part of 
the tensor within square brackets and 

(L/% = 2 V (i pj) - | gi j V fc (3 k . (2.16) 

Since we are focusing our attention on the method- 
ologies followed by CBB and GGB, we introduce at this 
point the assumptions that both share, namely conformal 
flatness g^ = rjij and vanishing of K. These two condi- 
tions exhaust six of the eight freely specifiable quantities 
in {gij, K^). Under these assumptions, Eqs. 1(2.9(1 and 
lUTfgll reduce to: 

V 2 </>+ i0 5 iyi y ' =0 (2.17) 
8 

Vjl y + 10 A ij Vj \n.(j> = . (2.18) 

In the remainder of the paper, we will continue refer- 
ring to the approach followed by CBB as CTT and by 
GGB as CTS, but one should keep in mind that CTT 
and CTS are more general. They do not necessarily re- 
quire conformal flatness nor K = 0. 

A. Conformal- Transverse- Traceless Decomposition 

The CTT method was pioneered by York and collab- 
orators Until recently, CTT has been the preferred 
method to construct initial data in numerical relativity. 
Given Eqs. ((2.17|) and 1(2.18(1 . one introduces the follow- 
ing conformal transformation and transverse-longitudinal 
decomposition of the extrinsic curvature: 

A ij = cj)' w A tj 

= <j)~ w [O lj + (LIU)* 3 ] , (2.19) 

such that Vj Q i] = and 6 l 4 = 0. We label CTT quan- 
tities with hats. 

Here 8 y represents the transverse part of A l K Since 
U is trace and divergence free, it contains the two 
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remaining freely specifiable pieces of the initial data 
(gtj, Kij). The longitudinal part of is given by 
(hwyi . It is this part of A 13 ' that is fixed by the momen- 
tum constraint 12.18fl . Furthermore, it is not difficult to 
show that Q t3 and (hW) 13 satisfy the following orthogo- 
nality condition: 



Si 



Q ij {hW) ijv fqd 3 x = 0. 



(2.20) 



For simplicity CBB assumed that U = 0. Therefore, 
substitution of lf2~T9j) into lf2~T7l) and (|2~T%|) yields: 



+ - cj}- 7 (1LW) 2 = 



Vj(Lwy 3 = o 



(2.21) 
(2.22) 



where we have introduced the notation (LI-F) 2 = 
(LW) lJ (LW)ij. In summary, the initial data prob- 
lem under the CTT approach reduces to solving first 
Eq. H2.22fl for W l and then using this solution to solve 
for 4> from Eq. 112.2111 

Bowen and York [22| found solutions to the momentum 
constraint 12.2211 representing single black holes with lin- 
ear momentum P l or angular momentum J\ Specifically, 

W i = — -(7 P i + n* rij P j ) (2.23) 
4r J 



W i = -—e^ k J,n k 



(2.24) 



with n l the unit normal of constant r spheres in flat 
space; namely, n l = (x, y, z)/r = (d/dr) 1 and r = \\x l \\. 
In terms of the solutions (|2.23|1 and 12.24JI . the extrinsic 
curvature A 13 takes the form: 



A 13 

#3 



3 

2T 2 " 

6 

"3 " 



2p( i n j) - (rf 3 - n l n 3 )n k P h 



(i j)H 



Jk ni 



(2.25) 
(2.26) 



Given the solution with linear momentum (|2.25|l . so- 
lutions representing non-spinning binary black holes are 
obtained via the following superposition: 



A 13 = A l3 \C\, p£) + A l3 (C%, P 2 fc ) 



(2.27) 



where C\ with A 



1, 2 denotes the coordinate location 
of each of the black holes. The explicit expression for 

i* J '(Cl> p a) is S iven b y ESUl with r A = \ \x* ~ C\\\ and 
n\ = {x l — C 1 a )/ta- It is not difficult to show that, with 

A tj given by l|2.27|l . the total ADM linear and angular 
momentum (about the origin of the coordinate system) 
are P l = P{ + P 2 4 and J, = e ijk (C{ Pf + C{ PI}), respec- 
tively. 

CBB specialized to the case of equal mass black holes, 
in the center-of-momentum frame of the binary system 
and in quasi-circular orbit. Thus, P l = P[ = — and 
Pi C l = Pi {C\ — C\) — 0. In addition, asymptotic flat- 
ness requires that in solving the Hamiltonian constraint 



(|2.21() one imposes the following Robin boundary condi- 
tion as r — > oo: 



(1 



(2.28) 



The differences between the procedures followed by Cook 
@ and by Baumgarte and Baker p| entered in 
the approach used to handle the black hole singularities. 
Cook Rl used the conformal-imaging approach. Baum- 
garte 121 and Baker on the other hand, applied the 
puncture method developed by Brandt and Brugmann 
[l4|. Constructing sequences of initial data sets reduces 
then to a 3-dimensional parameter space consisting of 
C = ||C l ||, P ee ||P l || and a third parameter. For the 
conformal-imaging approach method the third parameter 
is related to the radius of the throat of the black holes 
and for the puncture method to their "bare" mass. 



B. Conformal-Thin-Sandwich Decomposition 

The fundamental difference between the CTT method 
used by GGB and the CTS method used by CBB lies in 
the decomposition of the extrinsic curvature. Instead of 
using H2.19fl . one introduces 



A ij = 



-4 A ij 
2a 



(L/3) 4 



(2.29) 



where bars will be used to denote CTS quantities. 

The first aspect to notice is that, strictly speaking, 
(L/3) y /2a and u 1J /2 a cannot be viewed respectively as 
the longitudinal and transverse parts of A tJ . That is, 
they do not satisfy the orthogonality condition l|2.20(l . In 
Ref. [2^, the measure used in 12.20fl is modified to show 
that it is possible to arrive to a orthogonal decomposition 
of A 13 of the form 



.1" = (-)'•' 



1 

2a 



(LP) 1 



(2.30) 



with VjS 4 - 7 = <d l i = 0. The orthogonality condition in 
this case reads 

/ e l] -!- {LB)ij 2a^/rjd 3 x = 

J St ^ 

f & ij (LP)y ^jd 3 x = 0. (2.31) 

The next step in the CTS approach is to realize that 
the freely specifiable nature of Uij can be exploited to 
find quasi-equilibrium configurations. From the evolu- 
tion equation H2.13|l one has that 



d t gij = (JLf3)ij -2a A tj = u t 



(2.32) 



For quasi-equilibrium configurations defined by an ap- 
proximate helical Killing vector, it is then natural to 



choose = 0. 



5 



With Uij — 0, the Hamiltonian and momentum con- 



straint equations (|2.17|) and 1)2. 18[) read 
^ 2 1 1 ^ r (L/3) 2 = 



r 32 a 2 
V^(L/3) ij - (L/3) ij 



Vj In (a i 



= 0. 



(2.33) 
(2.34) 



Because of the particular form used to decompose Aij in 
(|2.29(l . a recipe for fixing a is also needed. Since one of 
the assumptions used is K — 0, it is natural to require for 
S t to be a maximal slice, namely dt K — 0. Substitution 
of K = dt K = into the evolution equation (|2.14|l yields 



lapse arises from imposing anti-symmetry with respect 
to the isometry. Because A iJ = (L/?) 1 -' /2 a, the condition 
a\s = forces one to impose (L/3) y |5 = in order to 
have regularity of the extrinsic curvature. The problem 
encountered by GGB was that the conditions l|2.41|l and 
(|2.42(l are not sufficient to guarantee that (L/3) l:, |,s = 0. 
CGG designed a procedure to achieve regularity on the 
throats; however, Cook [24j showed that a consequence 
of this regularization is the introduction of constraint vi- 
olations. 



III. CTT-CTS MAPPING 



V 2 a 



4a 



(2.35) 



from which one can solve for a. It is convenient to rewrite 
(JOS) with the help of Eq. ifO^ as follows: 



7 



32 (a< 



(L/3) 5 



0. 



(2.36) 



Notice from (|2.12ll and (|2.15|) that choosing = 0, 
dtK — and K = does not necessarily imply that 
dt<p = nor dtAij = 0, which would be necessary in 
order to have an exact helical Killing vector. 

The boundary conditions needed to solve Eqs. H2.33fl . 
(|2.34() and 12.35f> are that asr-^oo, 



n l Via 


(1 


-a) 
r 


(2.37) 




(1 


-$) 
r 


(2.38) 


(? 




(— V 


(2.39) 



with the orbital angular velocity of the binary and 
(d/d(p o y the Killing vector of the flat metric for rotation 
in the ip direction. In the GGB set up, the black holes 
are represented by two sheets of the Misner-Lindquist 
manifold (see .1] for details). That is, there are two 
throats S\ and S2 of radii a\ and ai located at coor- 
dinates C\ and C\ connecting the two sheets. Boundary 
conditions at these throats are derived from isometries 
between the two sheets. These conditions are that at the 
surface S of each of the throats: 



a\s 



2 r 



(2.40) 
(2.41) 
(2.42) 

(2.43) 



Here r denotes the radial coordinate distance from the 
center of the corresponding throat and n l its unit normal. 
Also, = h l j (3 J with hij — r/ij — Tij the induced 
metric in S and = n» (3 l . Condition (|2.4U|) on the 



The objective now is to find a prescription to trans- 
late between CTT and CTS data. Consider initial data 
(jjij, K,j) C TT and (g lj7 K tj ) C TS computed by the CTT 
and CTS methods, respectively. Without loss of general- 
ity, because of the conformal flatness and K = assump- 
tions, one can set: 



9ij\cTS = e 4 9ij\cTT 
Aij\cTS = Aij\cTT + 



2-2, 



(3.1) 
(3.2) 



with traceless tensor and e a function. When re- 



written in terms of CTT and CTS quantities, these rela- 
tions take the following form: 



2aS- 



e 2 (LW) z] + ai 



(3.3) 
(3.4) 



with bars and hats denoting CTS and CTT quantities, 
respectively. 

Since by construction, 



Vj{LWy j ^0 and Vj 



2ad>- 



e and Oy must satisfy the following condition: 
V.a ij + (LM/) y V,-e 2 = 0. 



(3.5) 



Clearly, if identical data are to be produced by CTT and 
CTS methods, one would need to show that e = 1 and 
that the only solution to (|3.5|l is a.^ = 0. 

The main difficulty in relating CTT and CTS data is 
that one not only has to establish a connection between 
their corresponding equations but also their asymptotic 
conditions as well as how the singularities are handled. 
Since CTT binary black hole initial data are computa- 
tionally less demanding to obtain (they only require solv- 
ing Eq. H2.21|) l. a starting point in checking whether CTT 
and CTS data are the same is to verify whether the CTT 
data satisfy the CTS conditions dt<jij = and d t K = 0. 
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A. From CTT to CTS 

We will concentrate attention on the steps to bring 
CTT data into a CTS form. That is, the initial data 
(<7jj, Kij) are given by the CTT method, namely c/ij = 
4 ? fa and Kij = (f>- 2 {1LW)i 3 , with (LW)ij given by fFSfy 
or (|2.26|) and (f> computed from 



1 



v 2 = --<?r 7 (Lwo 2 



(3.6) 



To simplify notation, we have dropped the hats used to 
denote CTT quantities. 

Satisfying the CTS condition d t K = is trivial since 
this condition is essentially an equation for the lapse 
function. Following a procedure similar to that used 



in deriving 



with Aij = 4> (IW)y instead of 



Aij — (f> i (L{3) i j/2a, one obtains 
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V 2 (a4>) = -(a4>)<j)^(JLWY 



(3.7) 



Asymptotic flatness imposes on Eqs. <|3.6[1 and (|3.7|) 

the following conditions as r — > oo: 



(1- 



ri 1 Vj(a i 



[l-(atf>)] 



respectively. These conditions imply that 



a0 

which in turn yield 



= l + - + 0(r- 2 ) 
r 

l--+0(r- 2 ), 
r 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



1 — a/r 
1 + b/r 



= 1 - 



2r 



+ 0(r" 2 ) 
0{r- 2 



(3.12) 



Above, a and b are parameters related to the ADM [f| 
and Komar [l^ masses, respectively. In GGB, the con- 
dition for quasi-circular orbits is found by equating these 
two masses. 

Next is the condition dtgij — 0. Once again, a deriva- 
tion similar to what was used in obtaining (|2.32|l yields 



(L/3) ij =2a<f>- 6 QLVT) y ' 



(3.13) 



Solving H3.13fl a $ it stands is conceptually difficult. 
There are five equations for the three components of the 
shift vector. Thus, solutions are only possible in the par- 
ticular case when the tensor a<f>~ 6 (LW) U ' is purely lon- 
gitudinal. A system of equations for f3 l that is not over- 
determined is obtained by taking the divergence of l|3.13|l : 



Vj(L/3)« = 2 [lW) ij Vj(a<p- 



(3.14) 



The potential problem of computing f3 z from Eq. (|3.14() 
is that, in general, yields 



(L/3) ij = 2a<j>- 6 (lW) lj 



(3.15) 



with VjU tJ 



0. A non- vanishing Uij would imply dtgij ^ 
0. Therefore, the only possible way of reconciling CTT 
and CTS data is if the solution to Eq. (|3.14() yields u lj = 
0. In Ref. [Sc} j we showed that u 1 ^ w for binary black 
hole data with punctures at separations close to ISCO. 

Given the asymptotic behavior of <fi and a, as r — > 00, 
Eq. H3.13fl becomes 



(h(3) lJ = 2 (LWf 3 



(3.16) 



Therefore, (3 l — * 2 W l + B l , such that B l is a conformal 
Killing vector, i.e. (L_B) U = 0. For example, in the case 
of binary data in a co-rotating frame, 



/3* -» 2 w l + n — 



d_ 

depo 



(3.17) 



In summary, given the vector W 1 , the initial data is 
completely determined once (|3.6f) is solved. The solution 
to this equation does not require knowledge of a and (3 l . 
The role played by Eqs. (|3.7() and l|3.13|l is to provide co- 
ordinate conditions that recast the data into a CTS form. 
The lapse function is fixed by Eq. (|3.7|l , which is derived 
from the condition for E t to be a maximal slice. The 
shift vector is fixed, on the other hand, by Eq. I|3.13|l . 
which is needed to impose the quasi-equilibrium condi- 
tion dtQij = 0. 



B. Komar Mass and Generalized Smarr Formula 

We need now to address how to fix the orbital angu- 
lar velocity £1 in l|3.17|l . In GGB, fixes the boundary 
condition that determines the initial data. Our case is 
different; the initial data (0, W l ) are obtained following 
the CTT method, so it is not until one solves for the shift 
vector that the issue of specifying a frame of reference 
emerges. If one chooses to solve for the shift with the 
asymptotic condition (3 % = 2W\ the data (4>, W\ a, i ) 
should be viewed as computed in an inertial frame of 
reference. 

One can transform these data to a rotating frame Q 
by simply changing the shift vector to 



(3.18) 



This transformation assumes that the inertial and rotat- 
ing frame are instantaneously aligned. It also implies 
that the lapse function, spatial metric and extrinsic cur- 
vature are unchanged. 

But the question still remains about how to fix f2. One 
would like in principle that 

l»=(^) =aN^ + ^ (3.19) 
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is a Killing vector, with t the time coordinate associated 
with the rotating frame. If such a Killing vector exists, it 
implies that the physical metric and extrinsic curvature 
satisfy 



dtgv 



= 
= 0. 



(3.20) 
(3.21) 



As pointed out by GGB , this requirement is too strong 
since energy would have to be pumped continuously into 
the system to keep the binary at a fixed orbit, thus vio- 
lating asymptotic flatness. 

Instead, one can require that in the rotating frame the 
conformal metric and trace of the extrinsic curvature sat- 
isfy 



d t gij = 
d t K = 0. 



(3.22) 
(3.23) 



Eq. (|3.22|l is equivalent to the thin-sandwich condition 
Uij = 0. Similarly, Eq. (|3.23|) is nothing other than the 



condition for S( to be a maximal slice. 

Given (EH2l and as well as K = 

equations (|2.13|l and l|2.14|l yield 



0, the evolution 



(L/% -2aA ij 



(3.24) 
V'V.a-ai' 1 ^ = 0, (3.25) 



respectively. Substitution of (|3.24() into (|3.25() . together 
with the momentum constraint V jA VJ — 0, yields 



V\V l a-A, ] P)={). 



(3.26) 



It is important to emphasize that this equation does not 
require conformal flatness to hold. 

Integration of (|3.26|> over E t excluding the black hole 
singularities (sources) yields 



1 

47T 



(Via- A ij /3 j )dS i + 1 = 



where 



1 



l = f- f {Via-Aijp^dS*. 

77 A=l,2 "' Sa 



(3.27) 



(3.28) 



Here 5,4 are 2-spheres surrounding the black holes with 
dS z oriented toward the interior of the black hole. 

The integral in (|3.27|) involving the gradient of the 
lapse has a physical interpretation. Given a time-like 
Killing vector the total mass (Komar mass) in a 
space-like hypersurface is given by 



111 



1 



Mkmr = f - <t> Vf n M N v dS , 

4 7T 



(3.29) 



with n M a unit vector in E t normal to a 2-surface S en- 
closing the sources. is the time- like unit normal to S t . 
In particular, one could take S to be at spatial infinity 



and £ M — > a (d/dt ) t *, with (d/dt )^ the time-translation 
Killing vector of the flat metric. One can then show j2j| 
that 



M KMR = ^- <t V,<« <!S' 
Air 



Thus, Eq. (|3~27jl becomes 



M 



KMR 



— <b A ti dS l + I = 0. 

4 7T ' 



(3.30) 



(3.31) 



One must emphasize that here Mkmr is meant to rep- 
resent the result of the integral (|3.30|) . It does not imply 
the existence of a global time-like Killing vector in E t . 
That is, in general dt4> ^ and dtAij ^ 0. 

The second term in (|3.31|) has also a physical interpre- 
tation. If P i -> fl(d/dip ) 1 , 



2VLJ A 



DM 



1 

47T 



A i:i (3 j dS i 



(3.32) 



with J adm the total angular momentum in S t 0. 
Therefore, 



M. 



KMR 



2ilJ A 



DM 



0. 



(3.33) 



One can then use this equation to solve for f2. However, 
Eq. H3.33fl does not fix in the sense of providing a quasi- 
stationary situation. Given any solution (0, a, (3 l ) to the 
system of equations (l3.6H3.l4jl . one can always solve for 

from (|3.33(l . The value of VL obtained in this fashion 
yields coordinate conditions (a, such that Eqs. (|3.22|) 
and H3.23fl ar e satisfied. It does not, however, necessar- 
ily imply that the vector Z** in l|3.19|l is an approximate 
helical Killing vector. 

As pointed out by Cook [24|, perhaps the most valu- 
able aspect of the work by GGB is to fix Q, by requiring 
that the ADM and Komar masses are equal. With this 
condition, Eq. (|3.33|) becomes 



M 



ADM 



2nj A 



DM 



7 = 0. 



with 



M ADM = - — 

A TT 



VicbdS 1 . 



(3.34) 



(3.35) 



Eq. (|3.34|) is also known as the generalized Smarr formula 


The idea behind constructing data with M A dm = 
Mkmr is motivated by Beig's observation [2(| that for 
stationary space-times these masses are indeed equal. 
Therefore, specifying fi to fulfill Madm — Mkmr would 
hopefully yield a co-rotating frame for quasi-circular or- 
bits that minimizes a suitable norm of dt4> and dtA^ |20| . 

If the shift vector solution to Eq. I|3.14|l is such that a 
residual transverse-traceless part u%j = dt gij is present, 
it is still possible to derive an expression similar to l|3.34|l . 
In this case, one needs to modify I to be [2(j 



- E 

4tt ^ 



A=1.2 



Sa 



(Via- A ij /3 j )dS i (3.36) 



-L J A ij ~g^ Uij ^/ld 3 x. (3.37) 
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Even for non- vanishing Uij, the condition Madm — 
Mkmr can be applied. It will minimize the norm of 
d t Qn : in addition to those of d t <\> and d t Aij [2£j|. In 
Ref. [19j , this methodology was used to construct binary 
black hole sequences with punctures. The results are in 
close agreement with those by Baumgarte [l2j . 

It is important to mention that in the CTS initial 
data construction by GGB, Eq. (|3.34[) is only used as 
a consistency check. The GGB methodology is such that 
part of the outcome is the rotating frame itself. As dis- 
cussed before, the angular velocity f2 of the co-rotating 
frame enters in the boundary conditions needed to solve 
for the shift vector. That is, O is adjusted to yield 
Madm — Mkmr- In our approach, on the other hand, 
the initial data is given by the CTT method, with fl a 
derived quantity. To fix Madm = Mkmr, one adjusts 
instead [2(| the parameters a and b in 1|3.10|) and i|3.11[l . 

Finally, if instead of j3 l — » Q, (d/dip ) 1 , one has that 
ft -> V (d/dx o y, with {d/dx o y the Killing vector of the 
flat metric along the direction x z , the integral in the the 
second term in (|3.31|l has a different physical interpreta- 
tion. In this case, 

-!- / iy ft dS l = 2V Padm , (3.38) 

with Padm the total linear momentum in S t [5j]. The 
Smarr formula reads in this situation 



M KMR -2VP A 



DM 



1 = 0. 



(3.39) 



IV. CTS PUNCTURES 



We will now explicitly consider CTS data involving 
punctures |14| . One decomposes 4> and a as 



such that 



P 
1 



p 
i 



02 
T2 



hi 



r-2 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



respectively. Above, a and b are parameters and u and v 
regular functions. The parameter a and b can be related 
to the "bare" masses of the punctures [2(j ■ 

Sub stitution of l|4"l"1) and (IO) into Eqs. (|2~33|) and 
(|2~3S|) yields 



For asymptotic flatness, the following Robin boundary 
conditions are imposed far from the holes: 



(!-") 
r 

(1-v) 



(4.7) 
(4.8) 



As pointed out in Ref. |2l| , the problem with punctures 
in the CTS method is that a natural choice is 



lim a = 1 

r — >oc 

lim a — — 

r-*r A 



c A < 0. 



(4.9) 
(4.10) 



Therefore, the lapse function would necessarily have to 
vanish at some internal boundary. In order to have reg- 
ularity of u and v, one would have to have that (L/?) y 
vanishes at least as fast as a at that internal bound- 
ary. Imposing the vanishing of (L/3) y in general requires 
more freedom than is available. One could in principle 
set ca < 0. This choice, however, yields solutions far 
from being quasi-stationary. 

If the combination of punctures and CTS seems to 
be intrinsically troublesome, it seems then puzzling that 
CTS conformal-imaging and CTT puncture binary black 
hole initial data sets exhibit good agreement as the sepa- 
ration of the binary increases. It should then be possible 
to translate, using the procedure described in the pre- 
vious section, the CTT-puncture data into a CTS form 
in spite of the presence of an internal boundary where 
a = 0. To demonstrate that this is the case, we con- 
sider single black hole solutions. For large separations, 
not only does the solution approximate that of a single, 
spinning black hole far away from the binary system, but 
it also approaches the solution of a single black hole with 
linear momentum near each of the black holes. Specifi- 
cally, we will consider Bowen and York [23] single black 
hole solutions using punctures and show that it is possi- 
ble to bring those solutions into a CTS form. 

From H4.5fl . I|4.6|l and (|3.13(l . the equations to be solved 
are: 



V 2 u 



V 2 v 



1 



1 



(W 



- p 7 (hW) 2 
8 [pu + 1)' 

7 - iqV + 1) — (LIT) 2 
8(pw + l) 8 q [ ' 

(pu + l)< q 



(4.11) 
(4.12) 
(4.13) 



V 2 w 



V 2 v 



32 a- 



■(L/3) 2 



l(5) (L ^ )2 



(4.5) 
(4.6) 



with (LWyi given by (|2~2l|) or (|2~2^|) . Notice that be- 
cause (LW) lJ is linear in the momentum, we will only 
consider perturbative solutions |27l l28j | to second order 
in the momentum for u and v and linear order for f3 l . 
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Therefore, 



where 



V 2 u 



1 f 



(4.14) 



V 2 V = I ( f_ .*ff (W) 2 (4.15) 



; (f+i) 6 
(f-5) 



(f+i) 



(4.16) 



where over-bars here denote scaling with respect to the 
parameter a (i.e. f = r/a, b — b/a, P = P/a and J = 
J/a 2 ). 



A. Spinning Bowen-York Black Hole 

For a spinning Bowen-York black hole [2^|. the only 
non- vanishing component of W 8 and (LW)jj are 



3J 



(LW)rj, = — sm 

where J is the total angular momentum. Therefore 

12 J 2 



(4.17) 
(4.18) 



(W) 2 = ^sm 2 , 



[Po(M) - iM/*)] > ( 4 - 19 ) 



where /i = cosf? and p are Legendre polynomials. Sub- 
stitution of (|4~T3|) into Eqs. (|4.14II4.16[> yields 



V u u = -- 
V 2 v 



3 J 2 r 



2(7TTF [jPo(m)_P2(ai)] (4 ' 20) 

21 J 2 r (r - b) 



2 (r + 1) 
= Qj r 



2 (r-b) 



(r + 1)" 



[Po(M) - (4-21) 
(4.22) 



To simplify notation, we have dropped the over-bars, and 
rescaled [3 V as ft? a. 

Regular solutions to (|4.20() - (|4.22ll arc 

u(r,li) = J 2 {u (r)P (p)+u 2 (r)r 2 P 2 (v)} (4.23) 
v(r,fi) = J 2 [vo(r)P (^+v 2 (r)r 2 P 2 (^} (4.24) 

-2 3(6 + 3) 



fff{r) = J 



(r + 1) 3 2 (r+1) 4 

6(26 + 3) (6+1) (1-6) 
5 (r + 1) 5 (r + 1) 6 10 



,(4.25) 



40u (r) 

20w 2 (r) 
40w (r) 



1 



1 



1 



(r + 1) (r + 1) 2 (r + 1) 3 
4 2 



(r + 


I) 4 


(r + 


1)5 




1 






(r + 


If 






(3 6 


-4) 


(3 6 


-4) 


(r- 


H) 


+ -. 

(r ■+ 


-I) 2 


(3 6 


-4) 


(3 6 

+ V- 


+ 31) 



(4.26) 
(4.27) 



(18 6 + 32) 10(6+ 1) 



20w 2 (r) 



(r + 1) 5 (r + 1) 6 
r (6 -6) - 66+1 

(>Ti? ■ 



(4.28) 
(4.29) 



These solutions have the following asymptotic form 
when r — > oo: 



1 / J 2 
1+ r 1+ 40 



a 



I- 1 - 

r 



6+1 



J 2 (5-36) 
40 



J 



10 r 3 
J(l-6) 



10 



■2iy ¥ 



(4.30) 
(4.31) 

(4.32) 



Furthermore, integration constants for /3 ¥ have been cho- 
sen such that pv(0) = 0. From Eqs. (|3~5S|) and ll3~30l) . 

the ADM and Komar masses are respectively 



J 2 

M ADM = 2 + ^7j 



M 



KMR 



6+1 



J 2 (5-36) 
40 



(4.33) 
(4.34) 



On the other hand, it is not difficult to show that near 
the puncture the integral /, defined by (|3.28|) . yields 



/ = - [1 + 6+ J 2 (6 Mo (0) + Uo (0))] 
J 2 (56-3)' 



1 + 6 

Similarly from l (3~32) l. 

Q = 



40 



10 



(4.35) 



(4.36) 



which is consistent with our requirement that asymp- 
totically (3 V — * Q. Clearly, the computed values for 
Mkmr, I and and £1 satisfy the Smarr formula Mkmr ~ 
2 il J + I = 0. In addition, to satisfy the condition 
Madm — Mkmr,: one needs to set 6=1, which in turn 
implies that = 0. 
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B. Boosted Bowen-York Black Hole 



where 



The only non- vanishing components of W l and (hW)ij 
for a boosted Bowen-York black hole 221 are: 



W r 
W 6 

(hW) rr 
(hW) r9 



P 

-2— cos 9 

r 

7 P 

--z sin (9 (4.37) 

4 r 2 v ' 

3 — cos 9 

3 P a 
cost/ 

2 r 

3 

--p cose 

-~P cos 9 sin 2 9, (4.38) 



where the linear momentum is assumed to point along 
the positive z axis and have magnitude P. 

With the above expression, the term (LV+) 2 takes the 
form 

(W) 2 = |4l 5P oW+ 4P 2(M)]- (4-39) 
2 r* 



The equations to be solved for u and v in this case are: 



V 2 u 
V 2 w 



/,2 ^[5Po(^+4P 2 ( A1 )] (4.40) 



[5Po(m) + 4P 2 (m)],(4.41) 



16 (r 
21 P 2 r 3 (r - 6) 
16 (r + 1) 8 

and for the shift vector 



rdr^)-do^ = 9 P ( ; 4 | r i)7 5) cosfl (4.42) 

r 2 r /? fl + d e (3 r = -3 P . r ^ ~ ^ sin . (4.43) 

(r + 1)' 

In the equations above, we have introduced the same 
scaling with respect to the parameter a as in the previous 
section. 



Regular solutions to (|4.40|) - (|4.43|l are: 



u(r, fx) 
v(r, /i) 
/T(r, 9) 

P 9 (r, 9) 



P 2 [«o(r) P (m) + «a(r) P 2 (/i) ] (4.44) 
P 2 [t.o(r)Po(M)+«2(r)P 2 (M)] (4.45) 
PR(r) cos (9 (4.46) 



PS(r) - sin 0. 
r 



(4.47) 



32 



80 ru 2 (r) = 
+ 



32 



voir) 





1 


2 


(r 


+ 1) 


[r + iy 




1 


I 


(r 


+ 1) 4 


f 5 (r + l) 5 




15 


132 


(r 


+ d + 


(r + 1) 2 




53 


96 


(r 


+ 1) 3 


f (r + 1) 4 




82 


84 


(r 


+ l) 5 


r (r + l) 5 


84 


ln(r 4 


■1) 


(6 




(6+15) 



(4.48) 



(4.49) 



(r + 1) (r + 1) 2 
(6 6 + 20) (8 6 + 15) 
(r + 1) 3 + (r + 1) 4 
(23 6 + 30) (6+ 1) 



80 r 2 v 2 (r) 



5 (r + 1) 5 (r + 1) 6 
(2016+ 1020) (275 6 + 825) 



(4.50) 



(r + 1) 



(287 6 



(r + 1) 2 
(185 6 + 283) 



(r + l) 3 (r + l) 4 

(66 6 + 80) 10(6+ 1) 



(r + 1) £ 



(r + 1) 6 



84 



+ — (6 + 8) ln(r + 1) - 105 



(4.51) 



R(r) 



1 



(r + l) 6 
1 



" 1 (&_5) + |(6-5)r 



(1516-751); 



80 

^(103 6- 503) r 3 

^(116-51)r 4 -4r 5 

lib-5) 



(4.52) 



Sir) 



1 



(r + 1) 6 
1 



I(6-5)-|(6-5)r 



(149 6- 749); 



80 

1(97 6- 497) r 3 
A(96- 49)r 4 + Ir 5 



+ o( & "5). 



(4.53) 
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As 



oo, these solutions have the following form: 



1 + 

1 - 
-P 



1 

6+1 

(6-5) 
8 



5P 2 



32 



5P 2 (7-6) 



32 

cos 9 



_ p (&5) COS + 2 vr 



(4.54) 
(4.55) 

(4.56) 



(6-5) 
8 



_7_ 

27 



sin# 



(b — 5) sin (9 „ 
P y - ' + 2W e . 



From 



8 r 

and l|4.57|l . one has that asymptotically 
(5-6) 



(4.57) 



■ P l + 2 W l 



(4.58) 



where W l is given by l|2.23|l . 

From Eqs. (|53oTl and (jOUjl . the ADM and Komar 
masses are respectively 



Madm 



5P 2 



16 



MfiTMfl = 6+1 
and from lpQ2^ . 



5P 2 (7-6) 
32 



K = ^P. 



At the puncture, it is not difficult to show that 
I = - [l + b + P 2 (bu (0) + v (0))] 



(4.59) 
(4.60) 

(4.61) 



1 



P 2 (36-5) 
32 



(4.62) 



It follows then that Mkmr, I and V above satisfy, as in 
the case of a spinning puncture, the Smarr relationship 
Mkmr — 2V P + I = 0. If in addition one imposes the 
condition Mad m = M KM r, then 6=1- 20P 2 /32 + 
0(P 4 ) and V = P/2 + (9(P 3 ). 



C. Approximate Killing Vectors 

In the previous two sections, we showed that it is pos- 
sible, at least at the perturbative level, to bring single, 
spinning or boosted CTT-punctures into a CTS form. By 
construction, the data satisfy the conditions dtgtj = 
and dtK — 0. Furthermore, in these cases the problems 
associated with the vanishing of the lapse are avoided. 
What remains is to investigate how close the data come 
to satisfying the conditions d t (j> = and d t Aij = 0, con- 
ditions necessary for the existence of an exact Killing 
vector. 








FIG. 1: Surface plot of Q = (dtA lJ d t A kl g lk 'g jl ) 1/2 for a spin- 
ning puncture in units of J 2 . 








FIG. 2: Surface plot of Q = {d t A l] d t A k ig lk g 3l f /2 for 
boosted puncture in units of P 2 . 



From Eqs. I|2.12|l and (|2.15() . we have that 



d t Aij = CpAi, -2a A lk A k j 



[aRij-ViVja] 



TF 



(4.63) 



(4.64) 



where the condition K = has been used. Because the 
solutions in Sees. II V Al and HVBl are up to second order in 
the momentum (linear or angular) for (j) and a and first 
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order for /3 l , one only needs to focus attention to dt<j> and 
d t Aij up to second order in the momentum. 

From Sec. IIV Al we have that for a spinning puncture 
the solutions are independent of the coordinate tp and the 
only non- vanishing component of the shift vector is 
It is then clear from Ij4.63|l that dt<j> = 0. For the case of 
a boosted puncture, on the other hand, one has that 

P If f>AC H 

dt<p= wVTW {b ~ m0r2+5r+1)+O{p3y (465) 

At the end of Sec. IIVBI we saw that in order to have 
Madm = Mkmr, one needs b = 1 + 0(P 2 ); therefore, 
dt<p = 0(P 3 ) w 0. To seoncd order in the momentum, 
dt4> vanishes identically. 

To zero order in P or J, Eq. 14.6411 equation reads 

dt = r^2 dia g[~ 2 ' 7-2 > 7-2 sin02 l - ( 4 - 66 ) 

r (r + \y 

Therefore, dtA^ = since to this order Madm = 
M KM r yields 6=1. To first order in P and J, Eq. 
vanishes identically. 

On the other hand, to second order in J, we have that 
for a spinning puncture, even for 6=1, the only van- 
ishing components are dtA rip and dtAg v . Similarly, to 
second order in P, these same components of d t Aij are 
the only ones vanishing, even if b = 1 — 20P 2 /32 as re- 
quired by Mkmr = Madm- The non- vanishing of dtAij 
is more apparent from Fig. ^ and [21 where we plot the 
quantity Q = {dtAijdtAkig^g^ 1 ) 1 ^ 2 as a function of the 
coordinates r and 9 for both cases. 

Having dt Aij ^ to second order on P and J implies 
that the Bowen-York puncture solutions do not have an 



exact Killing vector to this order. This is consistent with 
the results in Refs. [27l |2S| showing that Bowen-York 
initial data for a spinning or boosted black hole do not 
represent a constant-time slice of a Kerr or a boosted 
Schwarzschild black hole, respectively. 



V. CONCLUSIONS 

We have considered single boosted or spinning black 
hole punctures of the Bowen-York type and demonstrated 
that, at the perturbative level, it is possible to establish 
a direct connection between CTT and CTS initial data 
satisfying the conditions used by GGB needed for the ex- 
istence of an approximate Killing vector. Since for widely 
separated binary systems the black holes can be viewed 
as nearly isolated, our results contribute to explaining 
the agreement between CBB and GGB binary data sets 
as the separation of the binary system increases. 
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